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We study the behavior of two planes of Quantum Heisenberg Antiferromagnet in the regime 
in which a Chiral Spin Liquid is stabilized in each plane. The planes are coupled by an exchange 
interaction of strength J3. We show that in the regime of small J3 (for both ferromagnetic and 
antiferromagnetic coupling), the system dynamically selects an antiferromagnetic ordering of the 
ground state chiralities of the planes. For the case of an antiferromagnetic interaction between the 
planes, we find that, at some critical value J3 of the inter-layer coupling, there is a phase transition 
to a valence-bond state on the interlayer links. We derive an effective Landau- Ginzburg theory for 
■ this phase transition. It contains two 17(1) gauge fields coupled to the order parameter field. We 

I study the low energy spectrum of each phase. In the condensed phase an "anti-Higgs- Anderson" 

mechanism occurs. It effectively restores time-reversal invariance by rendering massless one of the 
gauge fields while the other field locks the chiral degrees of freedom locally. There is no phase 
^ , transition for ferromagnetic couplings. 
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I. INTRODUCTION 
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The discovery of superconductivity at high temperatures in the otherwise insulating copper oxides has motivated 
a thorough search for new physical mechanisms for both superconductivity and antiferromagnetism. This search has 
produced a host of new possible mechanisms many of which are not yet established on solid ground. Among these 
new ideas, the anyon mechanismtl stands as, perhaps, the most novel of them. For this reason, it has attracted a 
lot of attention. At a microscopic level, the anyon state requires that the underlying insulating state, known as the 
Chiral Spin LiquicB (CSL), should necessarily break Time Reversal (T) iavariance and Parity (P). An experimental 
signature of a state with broken T and P invariance is optical dichroismB. So far, however, there is no experimental 
evidence in support of the spontaneous breaking of either T or P in the copper oxidesa. Clearly, the simplest option 
is that these symmetries are not broken in the copper oxides and that the insulating states are unrelated to the CSL. 
£h \ At the present time this appears to be the case. 

In this paper we will explore the possibility that T and P may be broken in one individual plane but not on the 
system as a whole. Individual isolated planes may still be in states which break T and P but the sign of this breaking 
may not be the same from plane to plane. The simplest case is to imagine that the copper oxide planes are coupled 
by some interaction and that this .coupling is responsible for the selection of the state. A version of this problem has 
been studied by Rojo and LeggettO. They considered two planes with a doped CSL on each plane and, hence, had an 
anyon superconductor on each plane. They further assumed that the planes were coupled together only by a direct 
Coulomb interaction between the anyons on each plane. They did not fix a priori the relative sign of the statistics 
of the anyons on each plane but, instead, asked which relative sign was preferred by the Coulomb interactions. They 
found that the Coulomb interactions prefer the relative statistics to be antiferromagnetic ordered, namely opposite 
signs. The Rojo-Leggett result is due to a rather subtle edge effect. In fact, they found no effect in the bulk. 

In many copper oxides, the physical situation is such that the planes come in groups in which the planes are 
closer together than among nearby groups. This is rather common in the Bismuth based copper oxides. Because 
in these materials the inter-layer exchange constant which couples the copper spins can be comparable to the intra- 
layer exchange constant, there is a competition between intra and inter layer types of ordering. Quite generally, one 
expects to find to distinct regimes in the phase diagram for bilayers. At weak interlayer coupling, the ground state of 
the individual layers may be stable. However, if the interlayer exchange coupling dominates, the likely ground state 
should be a valence,bond state on the interlayer links. The case of two coupled Neel states was considered recently 
by Uhbens and LeeQ, by Millis and MonienEI and by Sandvik and ScalapinoB. These authors considered the effects of 
an inter-layer exchange interaction on the Neel ground states of the planes. 

In this paper we will reconsider the problem of a bilayer of quantum antiferromagnets in a regime in which there 
is enough frustration to drive each plane separately into a Chiral Spin Liquid. The planes will be assumed to be 
coupled by an antiferromagnetic exchange interaction of strength J3. The problems that we want to address are: 
(a) does the inter-layer exchange interaction select the relative ordering of the chiralities and (b) what is the phase 
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diagram for this system as a function of the inter-layer interaction. We consider a situation in which there is a CSL 
ground state on each plane, with fixed chirality but arbitrary sign. We find that quantum fluctuations around this 
state select an antiferromagnetic ordering of the chiralities. This is a rather interesting result. It means that even if 
on each plane the system was allowed to break P and T, the dynamics selects the state which is on the whole P and 
T invariant. We also find that, as J3 increases, there is phase transition to a state that we identify as a valence bond 
state on the inter-layer links, nameJy a T and P invariant spin gap state very similar to the one found by Ubbens et 
alB, Millis et au and Sandvik et aa. The problem of the ordering of chiralities by an inter-layer exchange interaction 
was considered previously by Gaitonde et au. By means of a perturbative expansion in powers of J3 they concluded 
that the chiralities order ferromomietiailly. The results that we report here disagree with those of Gaitonde et al. 

As it is by now well knowriEjEij&Ej, the CSL state and its low-lying excitations can be described in terms of 
an effective continuum field theory which is very much analogous to a set of Dirac self-interacting fermions in two 
space and one time dimensions. We find that the essential physics of this system can be understood in terms of the 
properties of an effective continuum theory of Dirac fermions on each plane provided that a physically sensible cutoff 
is introduced. The effective model contains two sets of massive Dirac fermions on each plane. The chirality of the state 
is given by the sign of the mass term. As in WWZ, the fluctuations around the CSL of each plane are represented by 
gauge fields (one for each plane) . By a detailed microscopic analysis we find that the interlayer exchange fluctuations 
are represented by a complex order parameter field. The effective theory is controlled by three parameters: I) the 
magnitude of the fermion mass on each plane (i.e., the fermion gap in the CSL), 2) the interlayer exchange constant 
(which determines the energy gap for fluctuations of the order parameter) and 3) the number of fermionic species 
(which we take to be N). In this picture, the phase transition to the valence-bond state becomes the phase transition 
to a state in which the complex order parameter acquires a non-vanishing expectation value. Our basic strategy is to 
first derive this effective theory and then use it to address the issues of the ordering of chiralities and of the nature of 
the phase diagram. 

Mean-field theories of frustrated antiferromagnets on a single plane have yielded a host of possible non-magnetic 
variational ground states. The actual phase diagram is not known in detail although it is generally accepted that 
non-chiral states are somewhat favored by variational calculations. In this paper we will not consider how interlayer 
couplings may alter this competition among possible single layer variational states. Rather, we will describe how 
interlayer interactions disrupt the CSL in favor of an interlayer valence bond state, which is clearly favored at strong 
coupling. The determination of the global phase diagram for bilayers is an interesting problem which is however still 
outside the reach of present theoretical tools and beyond the scope of this article. 

The effective field theory of fermions can be studied within a 1/N expansion. We use this expansion for two different 
purposes. First we look at the quantum corrections to the ground state energy of a system in which the two CSL 
are decoupled. We find that, at leading order in the 1/N expansion, the state with antiferromagnetic (opposite) 
chiralities is degenerate with the state with ferromagnetic chiralities. However, we find that the leading corrections, 
due to fluctuations of interlayer exchange processes, the state with antiferromagnetic ordering of chiralities is selected. 
In addition to the spontaneous breaking of this discrete symmetry (the relative chirality), the fermionic theory for 
the bilayers undergoes a dynamical breaking of the interlayer (out-of-phase) gauge symmetry at a critical value of 
the interlayer coupling constant. This phenomenon is strongly reminiscent of the breaking of chiral symmetry in 
the related (but not equal) field theoretic Gross-Neveu and Nambu-Jona Lasinio modelsEJ. Also, within this 1/N 
expansion, we find a phase transition from a regime in which the two planes have CSL ground states with opposite 
signs, to a state in which the inter-layer order parameter field condenses. We further investigate the physics of 
this phase transition by deriving an effective Landau-Ginzburg-type field theory, valid in the vicinity of the phase 
transition, i. e. for J3 ~ Jf . 

The degrees of freedom of the Landau theory, which is fully quantum mechanical, are the interlayer order parameter 
field and the gauge fields of the two planes. We present a qualitative study of the fluctuation spectrum of the two 
phases. The weak coupling phase has (almost) the same spectrum as that of two CSL with opposite chiralities: 
semions with opposite chiralities and gapped gauge fluctuations). However, a the phase with broken symmetry (in 
which the interlayer field condenses) displays an interesting "anti-Higgs-Anderson" mechanism: the condensation of 
the order parameter field causes a gauge fluctuation, which is massive in the unbroken phase due to the Chcrn-Simons 
terms, to become massless. This, in turn, implies that any excitation which couples to the gauge fields (the semions, 
in particular) to become confined by strong, long range, logarithmic interactions. The resulting spectrum of the 
condensed phase is equivalent to the low-lying spectrum of a ground state of local singlets, i. e. a valence-bond state 
on the interlayer links. The interlayer gauge field, remains massive and it effectively disappears from the spectrum. 
Thus, the "anti-Higgs-Anderson" mechanism wipes out all trace of broken time-reversal-invariance in the system. 
Unexpectedly, in this phase the system is actually more symmetric than in the non-condensed state. 

The paper is organized as follows. In section n we introduce the model for the bilayer and develop the mean field 
theory and briefly discuss the phase diagram. In section III we address the problem of the dynamical selection of 
chiralities. In section IV we derive a gradient expansion for the low energy modes of the (two) gauge fields and the 
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relevant (scalar) channel of the field coupling the planes. In section [v| we discuss the properties of the symmetric 
phase where the field coupling the planes does not condense, and an effecti ve a ction for the gauge fields is derived 
and studied. Section VI deals with the broken symmetry phase. Section VII is devoted to the conclusions. We 



also include appendices which contain technical details of the mapping onto the effective continuum theory and the 
computation of Feynman diagrams relevant for the phase transition, the ordering of the chiralities and the gradient 
and it expansions. 



II. MEAN FIELD THEORY FOR TWO COUPLED CHIRAL SPIN STATES 



Our model consists of two square-lattice spin- A Heisenberg antiferromagnets coupled through an exchange interac- 
tion of nearest-neighbors spins between planes with strength J3, and nearest-neighbors ( Ji) and next-nearest-neighbors 
(J2) interactions on each plane. The lattice Hamiltonian reads 

H = H L + H u + J 3 J2s L (S)-S u (x + e z ) (1) 

X 

where Hl,u is the usual Heisenberg Hamiltonian, 



Hl,u = Ji E S l,u{x) ■ S L}U (x + ej) + J 2 E Sl,u{x) ■ S L ,u(x + ei + je 2 ) (2) 

x,j— 1,2 x,j=-\- i — 

Using the slave fermion approach, the spin operator can be written in terms of fermionic creation and anihilation 
operators S(x) = cl c (x)<j a ^Cf3(x) with the usual constraint of single occupancy. We decouple the quartic terms by 
using a standard Hubbard-Stratonovich (H-S) transformation. Up to an integration over the H-S fields, the original 
theory is equivalent to the one that follows from the action given by the lagrangian: 

C = C L + C v - \- E IX»@)| 2 + E lcl(x)xMcu(x + e z ) + h.c.] (3) 

"3 

x x 

where 



C L = ^c* L (x) (idt + fi) c L (x) + ^p L {x) (c* L {x)c L (x) - 1) 

X X 

-j- E \^)\ 2 -t 2 E \x, L m 2 

x-j = l,2 x;j=+,- 

+E [ c l(x)X3,l{x)cl{x + ej) + h.c] +E [cl{x)Xo,l{x)cl{x + e x + je 2 ) + h.c] (4) 

x;j— 1,2 x:j— +,— 

where we have dropped the spin indices a, /3 to simplify the notation, with a similar definition for L\j. Here fi is 
the chemical potential and x means (r,t). The constraint of single occupancy is enforced by tha bosonic Lagrange 
multiplier field <p(x). This type of factorization was originally proposed by Affleck and Marstor£3 and by Kotliarliil. 
The H-S fields can be parametrized in terms of an amplitude Pj(x) and a phase Aj(x). This Lagrangian has a local 
symmetry if the Lagrange multiplier field tp transforms as the Aq component of a U(l) gauge field. 

The MFT consists in integrating out the fermions, at a fixed density, and treating the fields \j (x) within a saddle- 
point expansion. As it is well known, one serious problem with this mean field theory, is that there is-.no small 
parameter in powers of which to organize the semi-classical expansion. Following Affleck and MarstontJ, we will 
allow the number of spin species to run to N instead of 2, which is the case for the spin-i Heisenberg model. After 
re-scaling the coupling constant strengths J's and the fluctuating part of the fields, a 1-loop expansion of the fermionic 
determinant around the N — ► 00 Mean-Field solution can be performed by keeping the diagrams up to order j? . We 
have S e ff [ip, Xj] = NS [[<p, and the quantum partition function is Z = J T>xT>x*'Dpe lNS . 

There exists a whole family of solutions of the saddle-point equations. The simplest solutions are the valence bonds 
states and the flux phases. These may or may not be chiral. In this work we consider the problem of the selection of 
the relative chirality of a state in which there is a Chiral Spin Liquid on each plane. Thus, we choose a saddle point 
which represents Chiral Spin States on each plane and we will investigate which configuration of chiralities is chosen 
dynamically. 
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Wen, Wilczek and ZeeH (WWZ) have given a construction of the Chiral Spin State, which was first proposed by 
Kalmeyer and Laughlinu. WWZ begin with the flux phases, which have a uniform value for the amplitude of the n.n. 
H-S fields, say p{x) = p. This amplitude however, can fluctuate. The phases of the Bose fields on the n.n. links of an 
elementary plaquette have a circulation equal to ir or — n in mean-field value. This feature produces a collapse of the 
Fermi surface into four discrete points of the Brillouin zone: (±^-,±^-J at which two bands of states (positive and 
negative energy, "conduction" and "valence" bands) become degenerate. At these points, the excitation spectrum is 
linear and gapless. This allows for a mapping onto a discrete version of the Dirac theory with two massless fermion 
species of two-component spinors, with the "speed of light" equal to the Fermi velocity vf = 2ap. This gapless state 
can become unstable due to the effects of fluctuations. Several channels are known to be possible. If the staggered 
part of the fluctuations of the amplitude of the Bose fields on the n.n. links, picks up a non-zero expectation value, 
gaps will open up in the elementary excitation spectrum and they will provide masses (or gaps) to the Dirac-like 
fermionic excitations. These fluctuations can be seen to drive the flux phase into a dimer or Pcierls state and do not 
break time reversal invariance or parity. 

A mass term in a Dirac equation for a single two-component spinor Fermi field in 2 + 1-dimensions generally breaks 
T and P since the Hamiltonian, while hermitian, becomes complex. Since all three Pauli matrices are involved (two 
for the gradient terms and the third one for the mass term) there is no basis in which the Hamiltonian could be 
real. Therefore, the Hamiltonian is not self-conjugate and T is broken. However, in the case in which two species 
of fcrmions are present, the presence of such mass terms does not necessarily break P and T since they may have 
opposite sign for the different species. This is the case of the so-called Pcierls mass, which occurs in dimer phases. It 
is here where frustration comes to play a crucial role. By turning on n.n.n. interactions, WWZ allowed for additional 
H-S fields on the diagonals of the elementary plaquettes. The MF configuration for the phases can be arranged 
so that each triangle in an elementary plaquette is pierced by a flux equal to In this way, a time-reversal and 
parity breaking mass can be generated, i.e., one can provide a mass with the same sign to both fermion species 
in the plane. In order to perform the mapping onto the Dirac theory it is necessary to introduce four different 
field amplitudes at each unitary cell of four sites. This procedure can be done on the real space lattice by defining 
four sublattices and assigning an independent field amplitude to each one and expanding in gradients of the field 
amplitudestHI, or on the reciprocal latticeBLSI by expanding the lattice amplitude at each point as a linear combination 
of four independent fourier components amplitudes. On the reciprocal lattice these fields are the fourier components 
of the lattice amplitude centered at the four Fermi points. The low-energy physics of the system is determined by 
the scattering processes among these four amplitudes. Any of these procedures is equivalent to a folding of the first 
Brillouin Zone. 

In thejGSS, the mean-field ansatz for the amplitudes and phases of the H-S fields on the n.n. and n.n.n. links is 
given byEM 

Xi (e, e) = -xi (o, e) = Xi(e, o) = -\x (o, o) = ip 
e ) = -Xa(°> e ) = -X2i e > °) = x(°> o) = -ip 
X+(e,e) = x + (o,e) = -x + (e,o) = -~x + (o,o) = iX 

X- (e, e) = X- (o, e) = (e, o) = (o, o) = -iX (5) 

The fields Xj, with j — 1,2 or j = +, — are the H-S fields sitting on the n.n. and n.n.n. links respectively. The four 
different sublattices are denoted by (e, e), (o, e), (e, o), (o, o), where e and o mean even or odd site respectively. 

Once the mean-field H-S ansatz has been used into the Hamiltonian for one plane, convenient linear combination 
of the four field amplitudes can be arranged in the form of two two- component spinors and one can re- write the 
lagrangian for a single plane in the form of a lattice Dirac lagrangian with two massive fermion species. So far we did 
not include any fluctuations of the H-S fields. We will be interested in the fluctuating part of the phase of the H-S 
fields. 

In order to capture the physics of the system in the regime of long- wavelength, low-energy of the spectrum, we 
do not need the full lattice theory, but a linearized version around the Fermi points that keeps all the scattering 
processes that are responsible for the behavior of the low-energy excitations of the system. In the case of only one 
square lattice bearing a Chiral Spin State, . we .| aij rive to a 2 + 1-dimensional effective action involving two massive 
relativistic fermions coupled to a gauge fieldcni3ii3. The form of this action is given by 

S = J dx J dx 2 {tpi (i$ -4 - m i) ipi + $2 (0 -4 - m 2) ^2} (6) 

The continuum field tp a is related to the lattice amplitude ^> a by ipaix) = ^a(x)/a. We use a representation of Dirac 
gamma matrices in which 70 = 03; 71 = — ia-2 and 72 = —io~\, where o-j, j = 1,2, 3 are the usual Pauli matrices. The 
coupling to the gauge field (the statistical vector potential ) A M comes through the covariant derivative = — %$.. 
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The statistical vector potential is given by Aj = <fij/a = 2p(f>j/vF and Aq = p/vp, where <pj is the. fluctuating part of 
the phase of the Hubbard-Stratonovich fields on the n.n. links, ip is the Lagrange multiplier fielda'E3E3 and xq = vpt. 

The masses of the fermions come from the amplitude of the H-S field on the n.n.n. links and give a measure of 
the amount of frustration present in the system. These masses, although not necessarily equal in magnitude, have 
the same sign for both species. We assume that these amplitudes are fixed at their mean field values, since we are 
interested only in the effects of interlayer fluctuations. .— . I— . 

In what follows we adapt the methods of referencesB andtL 2 ! to the bilayer problemlla. We have a duplication of terms 
due to the inclusion of the second plane and new terms arising from the interplanar interaction. In the continuum 
limit, the action for the fermions in the low energy theory has two species of Dirac fermions on each plane coupled to 
both the intra-layer and inter-layer Hubbard-Stratonovich fields which mediate the interactions among the fermionic 
degrees of freedom. For simplicity we will assume the degree of chiral breaking is fixed and parametrized by two non- 
fluctuating masses mjj and mj,. These masses are given by m^^j = £Xl,u/ v f, being \l,u the mean-field amplitude 
of the Hubbard-Stratonovich fields on the n.n.n. links. We assume that the mean-field approximation amplitude of 
the H-S fields on the n.n. links p is the same for both planes. Consequently the Fermi velocity is also the same. The 
only low energy intra-layer bosonic degree of freedom left are the gauge fields of the upper and lower planes Ay and 
Al and the inter-layer fields Xz- 

The continuum action for the bilayer consists essentially of Eq. (|^) written twice with labels L and U for lower and 
upper plane and an inter-layer part given by the coupling between planes 

Smterlayer = J dx Q J dx 2 {ip L (<^o7ol + </?l7l T l + ^272^2 + ^3^) IpU + h.C.} 

-yj dx « J ^ 2 [[/(|^l 2 ) + ^(M 2 ) + ^(l^l 2 ) + ^(l^l 2 )] (7) 

In this expression ipL and ipu represent the two Dirac flavors f/'l't/ that live on the lower and upper plane of the 
bi-layer. The r-matrices mix Dirac flavors inside each plane. 

The intra-layer gauge fields, which represent intra-layer phase fluctuations on n. n. links have to be kept since 
they enter at the leading order in the continuum limit. There are other operators, with the form of fermion mass 
terms, that have not been included which do not contain any derivatives but they describe other types of intra-layer 
ordering which compete with the CSL. To include such effects would require a theory of the full phase diagram which 
is beyond the scope of this paper. 

The bosonicpart of the inter-layer action shown in the second line of Eq. (j?]) , comes from the corresponding bosonic 
terms in Eq. (g) 

S b = -j-Jdx J2 {\Xz{x)\ 2 + \Xz{S+e 1 )\ 2 + \ Xz {S+e 2 )\ 2 + \ Xz {x + e 1 +e 2 )\ 2 } (8) 

x— (e,e) 

where x is an even-even site on the lattice at, say, the lower plane. However, in going to the continuum limit it proves 
more convenient to introduce the rotation given by the linear combinations of the four H-S fields Xz(x) which link 
corresponding plaquettes of the planes 

(9) 
(10) 

(11) 
(12) 

In terms of the rotated fields, and after taking the continnum limit, the bosonic part of the action takes the form 

S b = — - [ dx 3 [tp* (x)(p (x) + (pl(x)(pi(x) + p* 2 {x)p 2 {x) + ip* 3 {x)ip 3 (x)} 
33 J 

= N J J&rf ^ Xip 'o(9) i Po(q) + \<p*(q)<pi(q) + \(p* 2 {q)^ 2 {q) + \<pt(q)<p 3 (q)} (13) 

In the second line of Eq. (y_3j) , Fourier transforms have been taken and the coupling constant g 3 has been re-scaled 
by in order to allow a^-expansion (see below). In other words A = y-, where <?' 3 = The fields ipy, 



cp (x) s 


d \ (Xz(x) 


+ Xz(x^ 


-ei) 


+ Xz(x- 


-e 2 ) 


+ Xzix + e i - 


he 2 )) 


<Pi(x) s 


* \ (Xz{x) 


~Xz(x^ 


- ei) 


+ Xz(x- 


-e 2 ) 


- Xz(x + ei - 


he 2 )) 


<Pi(x) s 


* \ (Xz{x) 


+ Xz(x^ 


- ei) 


- Xz(x- 


-e 2 ) 


- Xz(x + ei - 


he 2 )) 


<p 3 (x) s 


* \ (Xz{x) 


~Xz(x^ 


- e\) 


- Xz(x- 


-e 2 ) 


+ Xz(x + ei - 


^e 2 )), 
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j = 0, 1, 2, 3 also have been re-scalcd to p/vp. As a result, the effective coupling constant that controls the inter-layer 
fluctuations is 173 = 2a= = Jj,{2a) 2 jvp and has units of length. Throughout this work we use dimensions such that 
[h] = [e] = [vp] = 1 where h, e and vp are the Planck's constant, the unit of charge and the Fermi velocity 
respectively. We have a natural scale in our theory which is the lattice constant do, or the inverse lattice constant 
which we shall call A and characterizes the momentum cutoff. 

From the free part of the action, and the fact that we are working in 2 + 1 dimensions, it is clear that the dimension 
of the fermion operators must be A « (length) . The dimension of the operator tp is also that of A. The coupling 

constant g 3 is dimensional with [A] = ^ = A. This dimensional analysis tells us that the effective four-fermion 

operator which represents the interactions between the fermions of the two planes, is irrelevant at the weak coupling 
fixed point and that, if a phase transition exists, it should happen at some finite value of the inter-layer coupling. We 
will see that this is indeed the case. 

Now we integrate out the fermions and obtain the effective action 



>eff 



-iATrln 



tip L -m L ip 
p* - mjj 



+ S b 



(14) 



where we have defined 



f = </?3T3 + <A)70 + VlllTl + <P2l2T2 



(15) 



The saddle point equations are 



93 



-tf(0) 



dk 3 

(27T)'- 



-tr 



n M fc ) 
u Mo) 




(16) 



Formally, this integral diverges linearly with the momentum cutoff scale A. As in all theories of critical phenomena, 
we will absorb the singular dependence on the microscopic scale in a renormalization of the coupling constant. We 
can define a critical coupling constant g c as the value of the coupling constant at which the expectation values for the 
fields coupling the planes first become different from zero. Clearly the solution with < <p* >= is allowed for any 
finite value of the cutoff, no matter how large. This is the phase where the inter-plane field is not condensed. The 
non trivial solution will first occur at the value of the coupling constant g c - given by 



1 




tf-m L p(k) \ 
p*(k) t-mu ) 







Sip(k) 



MO) 




(17) 



evaluated at the point where the p's vanish. 

Notice that although the bare value of the coupling constants arc originally the same and equal to 33, they are 
associated with operators which do not scale in the same way. Their critical values are different as well. As an abuse of 
notation, from now on we are calling scalar to the interaction channel given by the field p 3 , frequency- vector channel 
to the field tpo and spatial- vector channels to the ones given by ipi and p2- 

Without coupling between the planes we have a degenerate situation between a state in which both planes have 
the same amount of frustration (i.e., the fermion masses are the same in magnitude) but their relative sign could be 
the same or opposite. We are going to call these two states ferromagnetic (FM) or antiferromagnetic (AFM) ordered 
respectively, understanding that we refer to the relative ordering of the sign of the chiralities. We want to investigate 
how the degeneracy between the FM and the AFM arrangement of masses is removed. For simplicity, we give the 
results for the case of \mp\ — \mu\ — m > 0. They may carry any sign. We define the variable s = sign(mi) sign(m[/), 
which takes values ±1. The critical values for the coupling constants are given by 



— = itr 

93 

1 

— = itv 
9o 

1 

— = itr 

9] 



S L (k)T 3 Su(k)r 3 
S i (fc) 7 oS c/ (fc) 7 o 
S L (k)'y j T j S u (k)'Y j 



^m(l-s) 

1 , 1 , 
40? 2tt v ' 



(18) 
(19) 
(20) 



where S a (fc) = -jj^- 



with a = L,U and j = 1, 2. 



6 



When the interaction between the planes is antiferromagnetic (i.e., J3 > ) the physical coupling constants remain 
positive. We are interested in the regime where m « A. For the case of an AFM relative ordering of chiralities (i.e., 
for s = — 1) we obtain 



1 

Hi 



A 



1 

7 



m 

7T 



1 



A 



rn 

7T 



(21) 



For m < A, we have c/ 3 < < g$] hence the channel which will first undergo a transition within the mean-field 
approximation, is the scalar channel, given by the field ^3. 

On the other hand, for the case of FM relative ordering of chiralities, we obtain 



1 A m — 1 A 

g\ 20? 7T g% g c j 40F 



(22) 



For m < A, again we have < g 3 < g? < g$. Again the channel which will first undergo a transition, if any, will be 
the scalar one. 

In the case of ferromagnetic inter-plane coupling (i.e., in the case J3 < 0) there is no transition, since the critical 
coupling constants always remain positive. The exact values of the critical coupling constants are not universal and 
they depend on the cutoff procedure that it is being used. Our continuum approximation is not very sensitive to 
these short distance features. However, the theory has a natural built in regulator since the model comes from a 
lattice theory. In other words, the qualitative feature of the existence of critical values for the coupling constants is 
independent of the type of cutoff procedure, although their precise value is not. The question of whether these critical 
values can be physically reachable is a different issue that needs a more detailed specification of the short distance 
properties of the model. We do not attempt to address this point here. We obtain the regularized saddle-point 
equations by subtracting the value of l/g c on both sides of Eq. (Jlq) . 
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(23) 



The simplest non trivial solution is the one where only the scalar channel ^3 is condensed. This channel has the lowest 
critical coupling, and it will be the first to pick a non-vanishing expectation value. For an antiferromagnetic relative 
ordering of the chiralities, which we will show it is favored in the case of antiferromagnetic Heisenberg exchange 
between the planes, we find 



7s " 7J ' 



dk 3 

(27T) 5 
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When solving Eq. (E4J) one gets 



(k^k^ - to 2 - |</? 3 | 2 ) 
yrr, 



k^ 



(24) 



(25) 



In Eq. (|25J), A 3 

is 



is the distance to the critical point. This is our equation of state. The non-trivial solution 



l^sl 2 =A 3 (A 3 -2m) (26) 

It is clear from Eq. ( p5|) that A3 < 0. When A3 < 0, i.e., when 53 > g 3 we find a phase where the scalar channel 
field has a nonvanishing expectation value given by Eq. (p6|). 

The physics of this state is the following. The fact that ip^ acquires an expectation value means that, on average, 
the inter-layer Hubbard-Stratonovich field is different from zero. Thus, it appears that in this state the fermions from 
one layer are free to go onto the other layer. However, the corrections to this mean field picture should, among other 
things, enforce the constraint of single occupancy at each site of each layer. The only state which is compatible with 
the single occupancy constraint and with inter-layer fermion hopping is a state in which, on each link between the two 
layers there is a spin singlet or valence bond state. Thus, the phase transition that we found is a transition between 
two CSL states on each layer (with antiferromagnetic ordering xif the relative chiralities) and a spin gap state with 
spin singlets on the inter-layer links. A number of recent worksQTj have predicted a similar phase transition in bilayers 
but between Neel states and spin gap state with properties which are virtually indistinguishable from ours. 
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III. RELATIVE ORDERING OF CHIRALITIES 



In this section we show that there exists a dynamical way in which the physical system selects a particular ordering 
of the chiralities in the planes. We assume that in each plane a CSS is stabilized. Thus, at each plane both Dirac 
fermion species are coupled to the mass term with the same sign. We assume that the mass is the same for both 
fermionic flavors in each particular plane, say mj, and mjj respectively. This is consistent with the fact that there is 
no explicit anisotropy present. As in section 0, the magnitudes of the masses are the same but their signs could be 
either the same or opposite. We neglect fluctuations of the n.n. amplitude of the H-S fields inside the planes, which 
can generate a difference between the masses of the Dirac species inside each plane, and even drive the CSS into a 
dimer phase (see, for example, referencelij) . 

Our goal is to compute the correction to the energy of the ground state of the bilayer system, due to the quantum 
fluctuations of the fields coupling the planes. We work in the phase where no field is condensed. Thus, the effective 
action derived in section |l| will describe the fluctuating part of these bosonic fields with zero expectation value. The 
strategy is, therefore, to expand this action in powers of (the small fluctuating part of) the fields (po to ip^ and keep up 
to the gaussian terms. Then, integrate the bosonic fields out and, after re-exponentiating the expression, obtain the 
desired correction to the ground-state energy density. This correction will contain a divergent part which is symmetric 
in the sign of the masses of the fermions in different planes, and a finite contribution which is a function of the fermion 
masses of both planes, rriL,u, with their signs. At this point, we look for the configuration of masses which minimizes 
the energy. The case of zero mass at any plane is excluded since we assumed beforehand that a CSS is stabilized 
at each plane. This is important since these masses provide the energy gap which is necessary for our saddle point 
approximation to be stable and to allow for a semi-classical expansion. 

The integration over the fermionic degrees of freedom gives the following contribution to the effective action (see 
Eq. (g)). We have 

= - i ™ 1 n(^— _ m °) + «f;^{(sQ)"} (27, 

v x n— 1 

Here 

At this point it is convenient to re-scale the fluctuating fields by ^7=. Under this transformation all the terms in 

Eq. ( p7| ) that are quadratic in the fields </?'s and Sb as expanded in Eq. ([l3]) become contributions of 0(1), being the 
classical energy of the ground-state (i.e., the classical part of the euclidean action) of O(N). To study the selection 
of the ordering of chiralities we need to compute this 0(1) correction to the ground-state energy due to the effect of 
the fluctuations of the fields coupling the planes. We first need to calculate the one-loop contribution to the fermion 
determinant. There is only one diagram to this order, which has two external bosonic legs and two internal fermion 
propagators 

i\ J dx 3 (sQSQ) = l 1 -J^J (s(k)Q( q )S(k - q )Q(- q )) 

As we saw before, this diagram has an ultraviolet divergence which will be absorbed in a renormalization of the 
coupling constants. So the kernels K^\q) in Eq. (p8| ) include both the finite part of the diagram and a contribution 
linearly divergent in the integration momentum. The computation of K, (q) , although rather cumbersome is fairly 
straightforward. Let us recall that we have four channels: <p3 can be regarded as a scalar-like coupling to the Dirac 
fermions; the other three -ipo to if2~ resemble a gauge-ficld-likc coupling. This is not the case, however, since Lorentz 
invariance is broken by the presence of the r-matrices in the expression for Q. This point is crucial. Since we do not 
have to preserve Lorentz invariance when regulating the divergent diagrams, time and space components do not enter 
on equal grounds. Our theory is in fact the continuum limit of a lattice theory. At that level it is very clear that the 
only physically sensible cutoff at hand is the inverse lattice spacing. As a result, our regulating procedure consists of 
integrating over frequency first and then using an isotropic gaussian cutoff for the spatial part of the momentum. In 
this way we expect to recover the qualitative features of the (finite) lattice theory in the continuum limit. Let us also 
mention that the only two spatially symmetric combinations of the interlayer amplitudes within a plaquette (see Eq. 
(|l2|)) are given by c/23 and tpo- 
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From now on, the expressions will be given in their Wick rotated (i.e., imaginary time) form. Consequently 
q 2 = q 2 + q\ + q%, where q Q — —iuu. We obtain (see appendix): 

/C(%)^-^{i(,W(l + S ))Z } (29) 



/C(0) ((?) = i_ {m(l + 2^ + 3^) +|- 



2 



2m 2 s 



4m 2 



3 4 
2 K ° 



- — 1 
2 V 



3Kn 



(30) 



/C«(g) = -4= - ^-{m (1 + 2k 2 + 34) + ( K . 2 - 4) 
4\/7r 2ir 1 y 3 2m v j 3 ' 



2 ° 



3 A q 2 



2m 2 s + Am 2 n 2 + ~k| - ^— (l + 3k£) 



(31) 



In Eq. ( J30| ) and Eq. (|3l| ) k 2 = q 2 /q 2 , with j = 0,1,2. Notice that the expression corresponding to the channel 
given by tpo (loosely speaking, the frequency channel) has an overall opposite sign to the expression for the channels 
given by Lp 1 and if2 for the finite part of the diagrams. However, the frequency channel does not have a divergent 
contribution. This sign will turn out to be quite important for the phase diagram. 
On the other hand, 2"o is 



Xq = — i sin 1 



v/4m 2 + q 



(32) 



At this point, in euclidean space, we have 



Z= / Vbe 



3 II / V<PjV<Pje*p\- 
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(27T) 



3 <p*(q) X-K.%(q) <pj(q) 



= #e- Eo exp ^ - ^ / In A - fC {j) (q) 

3=0 J 1 



(33) 



From Eq. (J33|), the correction to the energy of the ground state due to the fluctuations of the fields <p's is given by 
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2m z s + 4m z /eg + -«3 - ^— (1 + 3«5 



} 



d 3 g 
(2^)3 



In [A 



A 
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(34) 



We want to study the weak coupling regime, which corresponds to the case of large A in Eq. (|34|). Moreover, this is 
presumably the only regime for which Eq. (|3J) is valid, since as we show later, there is a critical value of the coupling 
constant at which there is an onset of condensation for some of the interaction channels between the planes. 
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By expanding in powers of j, to first order we obtain that the energy correction does not depend on the relative 
sign of the masses s and it is completely symmetric with respect to the exchange into mjj. This result remains 
true even when the magnitude of the masses are different. To second order we get 



AE^ = f. 



1 f d 3 q ( m 2 



X 2 J (2tt) 3 V4tt 2 



2 



4 

4m 



mq 



3 4m 2 + q 



2q 



1 



2b 



(35) 



The coefficient of s, where s is the relative sign of the masses (i.e., of the chiralities), is a function always positive. 
Thus, a minimum in the energy is obtained when s = —1, which indicates that the chiralities of the planes have 
opposite sign. This is the main result of this section. Recently, Gaitonde, Sajktar and RaoQ studied the problem of 
selection of the relative chirality by means of a perturbation theory in the inter-layer exchange coupling. They found 
that the ferromagnetic ordering was selected and that this result only appeared in third order in J3. This result 
disagrees with ours (see Eq. (p5|)). It is unclear to us what is the origin of this discrepancy. The work by Gaitonde et 
al relies on a rather complex lattice perturbation theory calculation of the inter-layer correlation effects. In our work 
we have evaluated the same correlation effects but within a continuum approximation which makes the computation 
more transparent and easy to check. We have used a cutoff only for the space components of the momentum transfers 
in our Feynman diagrams. The form of the cutoff that we chose closely mimics the effects of the lattice. Thus, it 
is unlikely that the discrepancy could be due to different choices of cutoffs. Similarly, the discrepancy appears at 
very weak inter-layer coupling where J3 <C \m\, where \m\ is the magnitude of the mass of the chiral excitations on 
each layer. Although it is conceivable that this discrepancy could be due to highly energetic processes which may 
be treated differently by both cutoff procedures, this appears to be unlikely since the mass \m\ is very large in this 
regime. Barring some numerical difficulty (which is possible in such involved calculations), the absence of a correction 
which depends on the relative sign of the mass in the lattice calculation (to the same order as the one given by Eq. 
([35|)) points to the occurrence of a special cancelation which we do not see in the effective continuum theory. We 
have also checked our result with other choices of cutoff on the space components and we have always found the same 
effect. Only in one instance, when we used a relativistic form of the cutoff, isotropic in both space and time, we found 
it necessary to go to third order in J3, which resembles the result reported by Gaitonde et al, but even in that case 
we found that the antiferromagnetic ordering of chiralities is the one energetically favored. However, the relativistic 
cutoff is certainly the one which is most unlike the lattice cutoff. In view of this considerations, we strongly believe 
that our treatment is robust and reliable. 



IV. LANDAU-GINZBURG EFFECTIVE THEORY 



We want to study the behavior of the low-energy modes for this system. The approach we are taking here is to 
derive an effective theory for the fluctuations of the cp-fields and the gauge fields. We want to study and characterize 
the phase diagram at the tree level approximation or Landau-Ginzburg approximation, and further on, investigate 
the effects of the fluctuations. We showed that there exist critical values for the coupling constants which possibly 
mark a transition between a symmetric or non-condensed phase for the </?-fields and a phase in which at least the 
scalar channel acquires an expectation value. The Landau-Ginzburg theory to be derived in this section will allow us 
to study the actual nature of this phase transition. We expand the fermionic determinant in a gradient expansion for 
slow varying modes of the fields in which we are interested. 

We derive an effective action only for the scalar channel. This particular channel is the one that first undergoes a 
condensation, for de case of an antiferromagnetic ordering of the chiralities, since it has the lowest critical coupling 
constant with a positive value. The other three channels will remain massive modes and consequently they can be 
integrated out of the theory. This process will involve renormalization of the parameters of the system but it will not 
affect dramatically the underlying physics. On the contrary, the scalar channel effectively undergoes a transition as 
the critical value of the coupling constant is approached and crossed. The bosonic excitations become massless at the 
transition point and we want to study the physics on both sides of this transition. We use the following definitions 
A^_(x) = A^(x) + A^j(x) and At(x) = A^(x) — Ay(x) for the in phase and out of phase gauge fields respectively. 
The covariant derivative is defined as T>^ = 8^ — 1 . 

The details of the calculation are described roughly in Appendix III. The following effective action is obtained by 
Fourier antitransforming the contributions of the 1-loop diagrams up to order -i, where N is the fermion species 
number. This includes bubble diagrams with up to four legs, since each of these legs represents the fluctuating part 
of either a matter or a gauge field, which has been previously re-scaled by a factor -^=. The loop integration adds a 
factor of N coming from the number of fermions propagating in the loop. From these diagrams we keep terms up to 
second order in the external momenta. In real space we find various terms; we get a contribution involving only the 
gauge fields which we call Sg^uge ■ This arises from the fermion loops corresponding to the propagation of spinon-hole 
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pairs inside each plane, without mixing. It contains the usual square of the field strength tensor and the induced 
Chern-Simons term. In the Ai_ ) coordinates this term is off diagonal, since the sign of time reversal invariance 

is opposite between the planes, 

(f^(x)F^(x) + F^F^ix)) (36) 



1 



647T 771 



dx 6 



The following term has a free part for the field ip and another part coupling this field to the gauge fields. A term 
coupling the gauge invariant current for the matter field <p to the field strength tensor of the in phase gauge field is 
also present. 

- afc W I ^^)^ J (-)W (37) 



In Eq. (37) we have defined the current operator for the field ip as 



j{-\x)= i {p*{x)d x p{x) - p(x)d x p*(x)) + JL A (-)( x)Mx) f (38) 

V TV 



Notice that all the terms are manifestly gauge invariant as it should be, since this symmetry was present before we 
integrated out the fermions. Notice also that the matter field couples only to the out of phase or relative gauge field. 
This is consistent with the symmetry of plane exchange which remains intact if the magnitude of the fcrmion mass 
is the same on both planes. In other words, the original theory was invariant under the exchange of Al and Ajj and 
the sign of the masses. This invariance should remain at this level for our approximation to be consistent. However, 
A/_) changes sign under this operation. This amounts to reverse the sign of the charge, or charge conjugation and 
consequently <p has to be conjugated. This renders the covariant derivative term and the gauge invariant current 
unchanged. On the other hand, F{*+\ is invariant under plane exchange. All the other terms are even on and 
our effective action verifies the plane interchange symmetry. Finally, from the contributions coming from the four leg 
diagrams which are of second order in the external momenta we can derive the following higher derivative terms 

S (2) (x) = Y^7r |mP^ J ^V*)! 2 - | fdx 3 ({VMp(x)T ({D»,Tr}<p(z)) 

We also get a self interacting term for tp given by 

S self (*) = - ^ ^ J ^\p(xr ~ Q - 1) / dx^(x)\ 2 (40) 

In Eq. (^p| ) above, we use the definition for g c introduced in section |l6|, i.e., -j- — In order to re- write this 

effective action in a simplified way we introduce some field re-scaling and define the following coupling costants, 

1 1 \ . 47r|m| 



\ [ dx a F^(x)F^JxMx)\ 2 + \ f dx 3 F(i\x)Fl i r ) (xMx)\ 2 } (39) 



cp(x) = y/Air\m\^(x); A^(x) = VN A^_y, m 2 Q = 4ar\m\ Q - iV A 



N 



N 2 167r|m| „ 1 1 
e = — G A = 7H — r5 and G = 



2tt' N ' 8\m\' 4m 2 ' 

By plugging all of these in, we obtain 

£ gauge = l(F^(x)A^ ) (x) + F^\x)A^ ) (x)) - + F 2 _ } (x)) (41) 

where Fx = e^xF^" is the dual of the field strength tensor. We also get a Lagrangian density for the field <j> given by 
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£ = d^*d^ - m 2 o \0\ 2 - A|0| 4 + d 



(42) 



where 



Cj = i {4>*{x)d^{x) - (f>(x)d^*(x)) A^ } (x) ~ G A F» +) {x) 



Af_ } (x) - 2G A F» +) (x)A^ (x)+g(1 F ( 2 _, (x) + \ 



(43) 



In Eq. (^) we dropped the higher derivative terms which appear in Eq. (|39|), except for the antisymmetric parts 
which involve renormalizations of the (different) effective charges for the in phase and out of phase gauge fields. 



V. SYMMETRIC PHASE 



In this section we want to study the physics of the regime in which there is no condensation of the field <jf> for the 
effective theory derived in the previous section, i.e., where <fi has a vanishing vacuum expectation value. This phase 
consists of the bilayer system with relatively opposite broken time reversal invariance between both planes, but the 
difference with the case of decoupled planes is that they are now linked through the fluctuations of the field <f>. This 
field represent a massive boson like mode with mass given by mg defined in the previous section by 



2 1 1 

rn = = 

9 9c 



Vp 
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(2a) 2 J 3 2y/n 2aJ 3 2^/tt 



(44) 



The magnitude of this mass measures the distance to the critical point. In this phase we are on the side of the 
transition in which mg > 0. It clearly corresponds to a weakly interplane coupling regime, i.e., the limit of small J3. 
The (/>-field can be integrated out to get an effective action for the gauge fields only. However, for our approximations 
to be consistent we need to assume that ?71q is much smaller that the fermion mass m. In other words, our results are 
valid on a window not to close to the phase transition (where </> becomes massless as g — > g c and m — * 0) but also 
not too far from the transition so that the mass of the collective mode represented by <j> never becomes comparable 
to the fermion mass. 

We are going to show that there is no renormalization of coefficient of the Chern-Simons terms that had been 
induced by the fermionic fluctuations on the planes, arising from the fluctuations of </>, at least to order . There are 
"charge" renormalizations in the sense that the coefficients of the field strength tensor for both ^4(+) and Ar_\ get 
renormalized. Furthermore, we will show that the spectrum of low energy excitations in this phase has two massive 
photons, whose masses do not violate the gauge symmetry but they break parity and time reversal invariance, and 
are very effective in taming the fluctuations of the gauge fields. In a sense we still have pretty much the same physical 
picture corresponding to two decoupled chiral spin liquid with opposite relative breaking of time reversal invariance. 
Consequently we will still have deconfincd spinon as the elementary excitations of the system. The issue of the statistic 
of the quasiparticles is a little more involved as we discuss below. 

Starting from the effective action derived in the previous section we can integrate out perturbatively the field <j>. 
This gives a result valid within the region of applicability of the gradient expansion. 

The integration over <j> gives the effective action (for small A^) We have, 



Vcj)V</>* expiS (4>,A^ 
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The cumulant coefficients can be computed in the usual wayEl to find 
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(46) 
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and 



(i {^d^ - <t>d x <P*)) (i - <R0*)) 



47T 



(47) 



Both integrals in Eq. dig) and Eq. ( J47| ) have a linear ultraviolet divergence and need to be regularized. One can 
use any of the usual regulators, for example Pauli-Villars or minimal subtraction (which is equivalent to an analytical 
continuation of the negative argument gamma function) However the finite part of both integrals after we treated them 
with the same regulating scheme is exactly the same but with opposite sign. This should be the case since it is required 
to preserve gauge invariance. In other words, we cannot generate a A^N 1 term in the symmetric phase because such 

a term would manifestly break gauge invariance and we know this is not the case. Therefore, the term A X A\ x ( \<f>\ 2 



in the r.h.s. of Eq. ( fi"5| ) should cancel exactly (and it does) the term A x A' 1 x (i 2 ((j>*d\<p — 4>d\4>*) (4>*d v 4> — cf>d v (f>*) V 
Notice that for the same token the term which could have given a renormalization of the cross Chcrn-Simons terms 
get canceled. In a sense, it is also gauge invariance which prevents the cross Chern-Simons terms to get renormalized. 

A minimal subtraction procedure will consist in the complete removal of the singular part. In fact, any cutoff 
procedure which preserves gauge invariance would work as well. It can be shown that our regularization prescription 
is entirely equivalent to the introduction of a gaussian spherical cutoff in the imaginary frequency (or euclidean) 
reciprocal phase space. A term of the form exp [— (-7-?){q% + m 2 )] does the job for us. One should be aware however, 
that this cutoff is not exactly the same used in sections || and [nj since there the cutoff was gaussian isotropic on 
the spatial components of the momentum but the frequency range was unbounded. Here that cutoff procedure would 
not work because it breaks gauge invariance. In sections [n] and III gauge invariance was not at stake and we were 
trying to implement a regularization that resembles closely what happens on a lattice. It should also be noticed that, 
although aparently the same field ip is involved in both cases, we were dealing before with ultraviolet divergences of a 
fermion loop integral, while here the field propagating is the bosonic field ip itself. In other words, we were dealing in 
the previous sections with the self-energy of the field tp while here we are dealing with the self-energy of the photon 
or the gauge fields. Finally we do get renormalizations for the F, 



(+) 



and F? 



terms. 



After this procedure is applied, we are left with the regularized (finite) form of Eq. (Eq) 
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(48) 



where 
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(49) 



In Eq. @ we used that F^ +) F x +) = 2F, 2 +) . 



Jjjq. < $io\) we useu niau (+) — (+)' 

We now explore the energy momentum dispersion relation. The low energy collective modes are fluctuations of the 
gauge fields. We will show that there exists a photon-like modeJsjut it is massive. This is of great importance for the 
survival of spinons in the energy spectrum (see for example ref.tHl) . The regularized (finite) theory has the form 



Seff (Ap) = / dx' 



- C-Fl_ } - c + F 2 +) + -e uvX (F^Al, + F^Al 



where c_ 



g-^v (+)(-) + (-)(+) 

64^1 m | (N - 2); and c+ = 64 ^ m | (N - §). In momentum space we have 



(50) 



Seff (Aa) 



A \ +) (-p) 



kqc^XvP J-i 



4 v 



= / A£(p) \c (p 2 g^ - P uPu)I a b + c 3 (p 2 g u „ - P uPu)T* » + 
Jp 

with a, b = (— ), (+) and /i, v the usual Lorentz indices, cq = — (c_ + c+), C3 = — (c_ — c+) and kq = —i 



(51) 
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This is a bilinear form in A% and the propagator for the gauge fields is just the inverse of the matrix shown 
in Eq. ([5l|). However, this matrix is singular unless we fix a gauge. This is so because the gauge field propagator 

I^A 1 ^ is not a gauge invariant operator and does not have a physical meaning unless we are working in a particular 
gauge. We need to add gauge fixing terms in order get the propagator. We may add for example, — i (dpA^_-^ 

and — i ^d^A^A , which in momentum space take the simple form —^p^Vv an d ~^PnPu- The three operators 

P^ u = p l _ l p Vl G^u = p 2 g^u and K^ u = e^XvP satisfy a closed algebra, and now the matrix can be inverted. After 
some lengthy though fairly straightforward algebra one gets 



n(-X-) n(-)(+) 

n(+)(-) „(+)(+) 



(52) 



given by 
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(55) 



where we defined 



M 



e 2 



ph 



4 \m\ 



(56) 



64c_c + V (N-2)(N-l/6) V /(1-A )(1 „_i_ ) 

as the "photon" massEH. To leading order in i we can rotate back to the Al, Ajj coordinates to get (in th Lorentz 
gauge a = f3 = 0) 
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and 
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(58) 



To next order in jt corrections we find additional off diagonal symmetric mixing terms. 



VI. BROKEN SYMMETRY PHASE 



In this section we want to study the phase where the matter field <f> condenses. Let us assume that we went through 
the critical point into the phase where m 2 , in Eq. (|44| ) becomes negative. From Eq. (|||) we now have another possible 

solution with finite (^j, which actually minimizes the energy. This is the usual non-trivial solution for a double- well 

effective potential of a 4 theory. When m 2 , becomes negative, the solution cj) — now becomes a local maximum 



instead of a minimun. The value of the new local mini mun can be obtained by minimizing Eq. (42) to be 0q = — -jj-, 
where we are using the definitions given in section IV . If we plug in this constant value of <fio , Eq. (S) becomes 



A 2 



(x) 2G A F» +) (x)Ap{x) + G Q F 2 _ } (x) + 1 



*?+)(*) 



(59) 
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As in the symmetric case we have 



4e 2 



Ff +) (x) + F^x) 



Now our effective action for the gauge fields reads 

Sett {An) = / dx 3 

where the new coefficients are 



*8 



(+) 



(60) 



(61) 



2A 



iV 



57T TO 



' ) 



(62) 



where now TOq < 



1 

4e^ 



2V 



647t|to| 



(63) 



and we define kq to be 



ACf) = — 2i 



4e 2 



blG A 



N 



,G 



6 647t|to| 



1 TOr 



1 TOr 



1 + -- 



3 TO 2 



4 V Am 2 ) 8tt V 4 to 2 



1 TOr 



(64) 



(65) 



We still need a gauge fixing term for In this way we recover the structure of Eq. ( |5l"| ) with minor changes. The 

photon mass has changed to 



= \/""" : M 1 -|^l 



4 ItoI 1 



2 TOc 

3 TO 2 



(66) 



The propagator for the out of phase field still corresponds to a massive field, which also has a longitudinal 

component 



£>(-)(-) = L 

» v 2 c+ 



PfJ, Pu 
p2 



1 



P 2 -M* h 



PjJ, Pv 

<t>l 



(67) 



n(+)(+) = 



2 c 4 




M 2 



(68) 



n(-)(+) -£»(+)(-) =f _ M. 

" 4 V ' 4to 2 / 4c_c 4 



(69) 



Eqs. (p7|))(|6^) 1 (|69|) give the propagators of the gauge fields in the condensed phase. By assumption m\ is a small 
parameter, since our approximation is valid for the vicinity of the phase transition where TOq is small (in units of 
the fermion mass) measures the distance to the critical point. Notice that, in this phase, the expansion in powers of 
jf has become an expansion in powers of this new parameter. This is consistent with our approximation because a 
gradient expansion amounts to an expansion in powers of an inverse (large) length scale, which in our case is set by 
the fermion mass or, in other words, by the spinon gap of the decoupled system. On the other hand, the "photon" 
mass is fairly large in this phase. 
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We conclude this section with a qualitative description of the excitation spectrum in the Broken Symmetry Phase. 
Let us look first at the gauge excitations. In looking at Eq. (|6^) we notice the existence of a massless gauge mode in the 
spectrum. The existence of this massless mode implies that any excitation which couples to the component of the 
gauge field experiences an effective long range (logarithmic) force mediated by the massless mode. In particular, the 
spinon excitations of the individual planes (which are semions in the unbroken phase) become permanently confined by 
the massless gauge fields. Recall that the logarithmic force is actually replaced-twraiConfinig potential due to the strong 
fluctuations of the gauge fields dominated by monopole-like configurationsc3li§l!Jj. In a sense, this makes fractional 
statistics unobservable since the quasiparticles which were able to bear it are no longer present in the spectrum. This 
spectrum is consistent with the fact that the statistical parameter 9 is not well defined anymore in this phase (it is 
no longer a topological number) since, as can be seen in Eq. ( |6l| ) it is now modified by a term proportional to the 
magnitude of the order parameter <f> . Also, in the broken symmetry phase, the time reversal breaking mass coming 
from the coefficient of the Chern-Simons term is no longer effective in controling the fluctuations of this particular 
mode. Actually, in this phase, the Higgs mechanism that takes place conspires to give a mass to the gauge field 
breaking spontaneously its gauge symmetry, while leaving the in phase field j4(+) untouched. In some sense the 
breaking of the phase symmetry enables the in-phase gauge field to become massless. 

This phenomenon is in striking contrast with the conventional Higgs- Anderson mechanism in which a spontaneously 
broken symmetry renders a gauge field massive. The remaining out-of-phase component is massive and its mass is 
huge (see Eq. (pq)), i. e. of the order of the fermion mass. This huge mass supresses the fluctuations of the field 
and, in this manner, it restores the broken time reversal invariance that was present in the decoupled bi-layer 
system. In particular, this spectrum implies that the only allowed fluctuations of the bilayer system are such that the 
chiralities of the planes become rigidly locked locally. Only in-phase, long wavelength fluctuations of the chiralities 
are allowed. Since the two chiralities have opposite sign, we conclude that, in this phase, there is a local cancellation 
of the chiralities of the planes. Hence, chiral fluctuations are eliminated from the physical spectrum. Recall 113 that if 
a Chern-Simons term were to be present, the monopole configurations would be suppressed and fractional statistics 
would become observable. This is precisely what happen in the symmetric phase. 

The spectrum that results from our analysis of the phase with broken symmetry is strikingly similar tO[-the spectrum 
of the bilayer system in the singlet phase discussed by Sandvik and ScalapinoB and Millis and MonienQ. In fact, we 
believe that the two phases are the same phase and that the broken symmetry phase is a phase with spin singlets 
connecting the two layers. 

VII. CONCLUSIONS 

In this paper we have reconsidered the problem of the selection of the relative sign of the chiralities of two planes 
with Chiral Spin Liquid states coupled via an exchange interaction. We found that the exchange coupling selects the 
antiferromagnetic ordering of chiralities and, thus, that T and P are not broken in bilayers. This result holds for both 
signs of the inter-layer exchange constant J3. Hence, even if each plane has a net chirality, the bilayer system does 
not. Such a system will not give rise to any unusual optical activity in light scattering experiments. We determined 
the phase diagram of the bilayer system and found a phase transition to a valence bond (or spin gap) state. Our 
analysis reveals the presence of an unusual "anti-Higgs- Anderson" mechanism which is responsible for wiping out all 
trace of broken time reversal invariance in the valence-bond state. In a separate publication we will report on results 
on the quantum numbers of the excitations and on the form of the wave function for the bilayer system. 
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APPENDIX A: COMPUTATION OF THE BUBBLE DIAGRAMS 



In this section we go through the computation of the correction to the energy of the ground state. The expressions 
given by Eq. (p9|) to Eq. (|3l|) are obtained from a one-loop diagram. In momentum space, 



JC(q) 



dk 3 T r [(tf + m L ) 4 + mu)<p*(-q)] 



(2tt) 3 



— <t) — mt 



(Al) 



with the definitions of section [nj. In computing the expression given above the following identities involving 2 + 1- 
dimensional Dirac 7-matrices will be important 



7a 76 = 9ab + ie-abd c ; tr (7a 7b) 
tr (7a7o7c7d) 



2g ah \ tr (7 a 76 7 C ) = 2ie a bc 
2 (g a bgcd + gacigbc — gacgbd) 



(A2) 



The integral in Eq. (Al) needs to be regularized, i.e., we need to cutoff the unphysical ultraviolet divergence due 
to the integration over momentum k. There is a natural cutoff in the original theory, which is the lattice spacing. 
However, in going to the continuum approximation we encounter the usual field theory divergences. 

There are several methods for regulating this type of integrals. The important point is that they should preserve 
the physical symmetries involved in the problem. For the case of a gauge field there are well established procedures 
such as the Pauli-Villars or the dimensional regularization methods. It can be shown that they preserve transversality 
(i.e., gauge invariance). In the case of the ^-fields, we do not have such a symmetry to preserve. In fact, not even 
Lorentz invariance is preserved. We have a length scale given by the lattice spacing, which in turn provides the 
momentum cutoff A that was mentioned in Sec.( On the lattice there is no cutoff for the frequency integral. Thus, 
our regulating procedure is as fo llows. First we perform a subtraction at the level of the integrands. This is equivalent 
to write the kernels in Eq. (Al) in the following way 



(A3) 



In this expression only the first term has an ultraviolet divergence. The term between brackets is convergent and can 
be calculated using standard methods. In the computation of ICj(0) we integrate first over frequency without any 
cutoff. This integral is still finite. After this step, we introduce an isotropic gaussian cutoff for the space directions of 
order A sa — where ao is the lattice spacing. However, once the integration over frequency has been performed, the 
finite contribution of the divergent term K.j(0) is independent of the particular form of the cutoff being used. 

The computation of K.j(0) with j — 0,1,2,3 is nothing but the calculation of the critical coupling constants 
performed in section (O). We have 



/C 3 (0)=ttr S L (k)T 3 Su(k)T 3 



1 



20? 



A 



2tt 



m (1 



(A4a) 



£q(0) = iti 



Si(fc)7oS(7(fc)7o 



2tt 



m (1 



(A4b) 



1 a 1 n s 
=A - — m (1 - s) 



2vr 



K.j(0) = itr S L (k)'YjT j Su(k)jjT j 

where j = 1,2. To get JC 3 (q) — £3(0) we need to integrate 

dk 3 (k 2 + m 2 )q- (q- k) + to 2 (1 + s)q ■ (2k - q) 



(2ir) s 



(k 2 + m 2 ) 2 ((k-q) 2 +m 2 ) 



(A4c) 



(A5) 



where k 2 , q 2 , k and q refer to the imaginary frequency rotated form of the tri- vectors k^ and q^, and s — sgn(mi) 
sgn(m;7). 

For the frequency channel given by ipo we get after rotating to imaginary frequency 




(2k q +q - (k-q)j (k 2 +m 2 )-q- (2k- q) (2fc| + m 2 (l - s)) 



(k 2 + m 2 ) 2 ((k - q) 2 + to 2 ) 



(A6) 
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For the channels given by ipi and <f2 we obtain expressions similar to the one for </?o with ko exchanged by k\ or ki 
in each case. The kernel for the spatial channels also have an oppo site sign to /Co (q) — ICq(0). 
We may write down the denominators in Eq. (A5) and Eq. (A6) in the form 



1 



De [m 2 



du(l -U) A 2 dXe -Xr 2 e -X[m 2 +u(l-u)g 



(A7) 



In Eq. (A7) we performed the change of variables A = a + f3 and u = [3 /{a + /3), we defined I = k — uq, and once 
again, for simplicity, we restricted to the case |mi| = \mu\ = m. After integrating over k (or rather, I) and A and a 
simple change of variables, Eq. ( A5) becomes 



1 



2tt 



For the vector channels we get 



1 



li = 7T Wl + s)-— sin" 1 



\J 4m 2 + q 



2tt ; 
1 

2^ 



Z f = ^Z9n rn (s + 2k] + 3/$ + 



2mq 



Am 2 + q 2 Kjy 



9jj 



2//r« -! 4///- ( k 2 + -k* 



(1 - 3k 4 



m 2 (l + s) + ig 2 



\q\ 



\q\ 



(A8) 



\J Am 2 + q 



(A9) 



In Eq. ( |Ag| ) gjj only indicates that the channels given by tpx and <p2 give a contribution with sign opposite to the one 
given by <po- Also, we have defined Kj =_jk an< ^ c = 1 + ^V- The expressions given by Eq. (|29| ) to Eq. (|3l]) can now 



be obtained simply by combining Eq. (A4) with Eq. (A8) and Eq. ( A9) 



APPENDIX B: GRADIENT EXPANSION 



In order to obtain an effective theory valid for long wavelength excitations, in momentum space we only need an 
expansion to the few lowest orders in the external momenta of the diagrams, since each external momentum will 
generate a space derivative in the Fourier antitransformed expression. Here we pursue further the expansion indicated 
by Eq. (p7|). However, instead of computing the loop diagrams exactly as we did in order to calculate the correction 
to the ground state energy in section III , we expand up to second order in the external momenta. Notice that we are 
not integrating out neither the gauge fields nor the y-fields in this case. 



1. Diagrams with two external legs 



We showed in section [H that the classical energy is of 0{N). From the expansion of the logarithm of the determi- 



nant, to second order in power of the (small) fluctuating fields contained in the operator Q of section III which now 
also includes the gauge fields in the planes -or to 0(1) in the jj expansion-, we have 



i- I dx 3 
2 



SQSQ) =- 



tr 



qJk 



tr 



+ 2 x 



>q Jk 
i 



S L (k)i L (q)S L (k + qyL L (-q) 
Suik^uiqlSuik + qttui-q) 
tr S L (k)ip(q)Su(k + q)<p{-q) 



(Bl) 



As a shorthand we have used the notation J = J k . 

For the gauge fields alone, the diagram has an ultraviolet divergence that needs to be treated. We use dimensional 
regularization to ensure traHsversality, i. e. to preserve gauge invaria nce. The calculation is similar to the one shown 



(B2) 



in chapter VII of referenced. The first two lines of the r.h.s. of Eq. ( |Bl| ) give / U^(q)Al(q)A" L (-q) where 



H^(g) = * e^ X q X U L A L {q 2 ) + {q»q„ - q 2 g,„) Tl L s L (q 2 ) 



q ~~A» 
2\ 
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with 



and 



1 1 



2tt 



2|to| 



: sinh 1 




(B3) 



4m 



sinh 1 




(B4) 



and the corresponding expression (UU) for the upper plane. This is the full expression for the polarization tensor for 
the gauge fields. The small momentum limit for Eq. (|B2[) is given by 



An TO 



1 1 ( i \ 

r- T {q^ -q 9^) 

l07T to ' 



Notice that the fermion mass with its sign enters the antisymmetric part of TLt£. Therefore, for an antisymmetric 
ordering of chiralities in the ground state, U^ L and will bear opposite signs. In fact, the ratio t^t in our case is 

actually always positive 1, since the sign of the masses have already been taken into account when definig Sl and Sjj 
in section III, and we are considering the case when the magnitude of the masses is the same on both planes. From 



now on, we give the details for the computation of the 1-loop diagrams involving exclusively the scalar channel. The 
third line in Eq. (Bl) can be rewritten as 



S L (k) r 3 S v (k + q) t 3 | ip(q) | 



(B5) 



This expression is similar to the ones we encountered in the previous appendix. Since we are interested in a small 
external momentum expansion, the exact expression shown above can be approximated by 



2 i I / tr 

Jq Jk 

which gives the following result 

l<^(<?)| 2 



S L (fc) J2(-l) n (Su(k)A n Su(k) 



M?)| 5 



Air (to 2 )!/2^ q " 



(B6) 



(B7) 



The above expressions are valid up to second order in the external momenta of the 1-loop diagram, i.e., in the 
momenta of the field ip or in the momenta of the gauge fields. 



2. Diagrams with three external legs 



The next term in the expansion of the logarithm of the fermionic determinant is 



1 



N 



(sq)' 



3 i 
i= x - 

3 i 
-= x - 

N 3 



' q J k 

II 

f p J q J k 



Tr 
Tr 



S L (fc) 7M S L (fc +p)0(q)Su(k +p + q) (<p(p + q))*\ A» L {p) 
Suik^Suik + p) m-q))* S L (k+p + q)0(- P - q)} AUp) 



(B8) 



Both terms in the r.h.s. of Eq. (B8) are similar. I can easily be shown that the zeroth order term in external momenta 
vanishes. For the scalar channel only, up to second order in the external momenta, the three-leg 1-loop diagrams give 
the following contribution 
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-m-riLI {^-^){Ai[p)-A^p))5{p+ s + P )^ q )( v {-s)r 

v jV 4tt \m\ j gjSjP 

(A£(p) + A£(p)) p" s A 5 (p + a + p) <p(q) (<p(-s)Y 



1 1 

JV 87r |m| 5 



(B9) 



3. Diagrams with four external legs 



The fourth order in the expansion of the logarithm of the fermionic determinant gives 

4 



AN , 
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+ N 

i 

+ N 
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+ N 
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+ 2N 



d.r: 



(sq) 



Tr 



l,p,q,k 



Tr 



l,p,q,k 



l,p,q,k 



S L (k) ltl S L (k+p) lu S L (k +p + l)<p(q)8 v (k - s) (0(-s))*\ A£(p) A»(l) 
S u (k)^S u (k+p)^S u (k+p + l)if(q)S L (k - s) (0(s))*] A£(p) A»(l) 
S L (k)^S L (k+p)<p(q)Su(k+p+ q) lv Su{k - s) (£(-«))*] A£(p) A£(Z) 
Tr rs L (fc)£(p)Sc(fc + p) (#(g))* S L (fc +p + q)${l))Su{k - s) (£(-»))*" 



Tr 



terms involving four gauge fields 



(BIO) 



From this expression we are going to consider only the four first lines on the r.h.s. of Eq. ( BIO ) as they will show to 
be the relevant terms for our g radient expansion. As we did before, we consider only the scalar channel. The third 
term on the r.h.s. of Eq. ( B1C ) gives a total contribution, valid to first order in the external momenta that looks 



If T<b (^l 85 - + R tr ^ 7 ' ¥+b]yi(p)AUi)w{q)M-s)T- 



The contribution coming from the first term on the r.h.s. of Eq. (BIO) is 



The second term on the r.h.s. of Eq. ( B1C ) gives 



1 1 

iV T&7T (TO 2 ) 1 



(Bll) 



(B12) 



(B13) 



The origin of the relative sign between the antisymmetric parts of Eq. (B12) and Eq. (B13) is the relative sign of 
the fermion masses on the planes. 
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